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Abstract. We perform an analysis of preliminary data on hadron yields and fluctuations within the Statis- 
tical hadronization ansatz. We describe the theoretical disagreements between different statistical models 
currently on the market, and show how the simultaneous analysis of yields and fluctuations can be used 
to determine if one of them can be connected to underlying physics. We perform such an analysis on 
preliminary RHIC and SPS A-A data that includes particle yields, ratios and event by event fluctuations. 
We show that the equilibrium statistical model can not describe the K/n fluctuation measured at RHIC 
and SPS, unless an unrealistically small volume is assumed. Such small volume then makes it impossible to 
describe the total particle multiplicity. The non-equilibrium model, on the other hand, describes both the 
K/n fluctuation and yields acceptably due to the extra boost to the tt fluctuation provided by the high pion 
chemical potential. We show, however, that both models significantly over-estimate the p/vr fluctuation 
measured at the SPS, and speculate for the reason behind this. 
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1 Introduction 

One of the main objectives of heavy ion physics is to study 
the collective properties of QCD matter. It's equation of 
[ state, transport coefficients and phase structure, and the 
. dependence of these on energy and system size. 

Thus, the natural approach to study soft particle pro- 
, duction in heavy ion collisions is through statistical me- 
■ chanics techniques. Such an approach has a long and il- 
I lustrious history [Tl[51[3111] . A consensus has developed that 
the statistical hadronization model can indeed fit most or 
all particles for AGS,SPS and RHIC energies t5 , 6 ,7..8, 9. 

mm- 

The statistical model obtains particle yields by as- 
suming entropy to be maximized given the constraints 
imposed by energy and quantum number conservation. 
These constraints can either be imposed rigorously, as 
required for closed equilibrated systems, or on average, 
as required for a sub-system equilibrated with an unob- 
served "bath". Full energy and quantum number conser- 
vation is usually referred to as the micro-canonical ensem- 
ble, while the Canonical (C) and Grand-canonical (GC) 
ensembles assume that, respectively, energy and other con- 
served quantities can vary via system-bath exchange. In 
this work, we shall concentrate on the GC ensemble, as we 
see it as most appropriate for describing the statistically 
hadronizing fireballs produced in heavy ion collisions. 

Our approach is not universally agreed on by the heavy 
ion community; In fact, noteworthy attempts were made 
to explain the dependence of certain observables w.r.t. en- 
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ergy and system size (The "horn" structure in the K'^ /n^ 
ratio dependence on energy [12] , as well as the strange par- 
ticle enhancement at larger energies and system sizes |131 
[T^ ) on the limits reached in different regimes; Averages 
of observables over all events are the same in all three 
ensembles in the limit of either a "large" system (high en- 
ergy and system volume) or an observed small sub-system 
of a large system (mid-rapidity), away from these limits 
(low energy, p-p and p-A collisions), the three ensembles 
give very different results, and, given the low multiplicity 
of strange particles in these collisions, it is necessary to 
enforce exact conservation to correctly count the particle 
states (Canonical or micro-canonical ensemble). The onset 
of strangeness enhancement [15], and perhaps structures 
such as the horn, have been speculated to arise from the 
transition between these limits. 

Both the horn [16^ and strangeness enhancement [5] 
are also regarded as indicators of a deconfined state of 
matter. If this is correct, than the energy and system size 
where they develop indicates the point at which the sys- 
tem can be described as an equilibrated deconfined Quark- 
Gluon Plasma (QGP). 

Distinguishing between these two explanations is a phe- 
nomenological question: The ensemble assumption is a fal- 
sifiable one, in the sense that a physically inappropriate 
ensemble is very unlikely to describe the available experi- 
mental data, if the data includes yields and event-by-event 
fluctuations [T7l[T8lfT9] . Hence, a statistical model analysis 
at different energies and system sizes, of both yields and 
event- by-event fluctuation observables, will isolate the sta- 
tistical physics (if any) responsible for the observed fea- 
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tures. The present work is a step in this direction, focus- 
ing on a Grand Canonical analysis of top energy RHIC 
(v^ = 200 GeV/A Au-Au coUisions) and SPS (^/s = 17.6 
GeV/A Pb-Pb collisions) heavy ion data. 

Yields and fluctuations in the GC ensemble are calcu- 
lated via the quantum statistic "ideal gas" formulae 



(27r)3 x-^^V^+^^'^ ±1 



(1) 
(2) 



Where gi is the degeneracy, V the system volume, the up- 
per sign is for Fermions and the lower sign for bosons. It 
is sometimes confusing that states in a "strongly inter- 
acting" system can be accurately predicted via the ideal 
gas ansatz. The key insight [3] is that accounting for all 
strong excitations (resonances) is equivalent to counting 
the QCD "energy levels" . Since hadronic interactions are 
resonance dominated (a consequence of the confining na- 
ture of QCD) this is a good approximation. Lattice studies 
have confirmed the validity of this approach quantitatively 

The final state yield of particle i is computed by adding 
the direct yield and all resonance decay feed-downs. 



(3) 



all j — ■>% 



all j^i 

all j — >i 

The parameter A; corresponds to the particle fugacity, re- 
lated to the chemical potential by A^ = e'^'/-^. Provided 
the law of mass action holds, it should be given by the 
product of charge fugacities (flavor, isospin etc.). It can 
be parametrized it in terms of equilibrium fugacities A"** 
and phase space occupancies 7^. A hadron i with q{jj) light 
quarks, s(s) strange quarks and isospin 13 has fugacity 



^« Iq Is 



J3 



(5) 



The temperature and chemical potentials can be obtained 
from data by doing a fit [H] (The chemical potentials 
for strangeness. As, and isospin, A/3, are usually obtained 
by requiring strangeness and Q/B to be conserved). 

If the system is in chemical equilibrium then detailed 
balance requires that 7^ = 7^ = 1. Assuming 7g = 1 and 
fitting ratios gives the T 160 — 170 MeV freeze-out tem- 
perature typical of chemical equilibrium freeze-out models 
at SPS and RHIC piTllS]. Some fits [TPTT] also allow for a 
kinetically out of equilibrium strangeness quantum num- 
ber, which is found to be at chemical equilibrium (7^ ~ 1) 
at RHIC and slightly below equilibrium (7^ < 1) at SPS. 

In a system expanding and undergoing a phase tran- 
sition, however, the condition of chemical equilibrium no 
longer automatically holds, so one has to allow for the 



possibility that both 7s and 7^ 7^ 1. In particular, if the 
expanding system undergoes a fast phase transition from 
a QGP to a hadron gas, chemical non-equilibrium ^ and 
super-cooling [22] can arise due to entropy conservation: 
By dropping the hadronization temperature to 140 MeV 
and oversaturating the hadronic phase space above equi- 
librium (7q ~ 1.5,7s ~ 2), it is possible to match the 
entropy of a hadron gas with that of a system of nearly 
massless partons [5]. These are exactly the values found 
for 7g and T at SPS and RHIC in fits where 7^ was a fit 
parameter [^[TU]. 

The "7g = 1" and "7, fitted" approaches are different 
models of how hadrons are produced in heavy ion colli- 
sions. They are both based on statistical mechanics, yet 
differ in the physics underlying it. They are, in principle, 
distinguishable experimentally but have not been to date. 

Fig. 1 in [9] illustrates why an experimental test of 
either of these models is non-trivial: An increase in tem- 
perature acts in a very similar way on particles and anti- 
particles as an increase in 7: The abundances of both go 
up. The amount in which relative particle abundances go 
up is different, making the two models distinguishable by 
yields alone. However, this difference is not enough to con- 
vincingly disentangle the two models, given the large ex- 
perimental errors in measurements of yields and ratios. 

To falsify one of the models it is necessary to con- 
sider event-by-event fluctuations as well as average parti- 
cle yields 23 . If jq becomes as large as claimed in [9] , the 
pion chemical potential Att approaches e™"/"^. Near this 
limit (corresponding to the critical density for pion B-E 
condensation) it can be shown [23 that (A^tt) converges 
but ((ZIA^tt)^) diverges as 



{{ANr) 



-1/2 



(6) 



where e — 1 — Xt^c^™^^"^ . thus, while T and 7^ are cor- 
related in yields, they are anti-correlated in fluctuations. 
The measurement of a yield and fluctuation constrains 
both to a high precision. At least one of the models will 
be ruled out by inclusion of a fluctuation in the fit. Inclu- 
sion of more than one fluctuation can be used to test the 
remaining model. 

The question of chemical equilibration is is closely re- 
lated to the accounting for hadronization when modeling 
freeze-out. While a quantitative treatment of this topic is 
still lacking, it is obvious that if the freeze-out temperature 
is higher (i.e., 7^ ~ 1), the effect of hadronic interactions 
between chemical and thermal freeze-out is greater than if 
the system freezes out from a super-cooled state (7^ > 1). 

It is also unclear what, if any, is the effect on ob- 
servables of the evolution between hadronization (the for- 
mation of hadrons as effective degrees of freedom) and 
freeze-out (the moment when all hadrons decouple) |24j . 
As shown in [25j , the discrepancy with experiment would 
lessen if chemical and kinetic freeze-out were to coincide. 
While, as claimed in [25], such a high freeze-out tem- 
perature would spoil agreement with particle spectra, fits 
based on a single freeze-out model prescription have shown 
that, provided a correct treatment of resonances is main- 
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tained, particle spectra are compatible with simultaneous 
freeze-out p6 l[27l[28] . 

A way to distinguish between these scenarios is to di- 
rectly measure the abundance of resonances. Here, the sit- 
uation becomes even more ambiguous: As pointed out in 
[29] , resonance abundance generally depends on two quan- 
tities: to/T, where m is resonance mass and the chemical 
freeze-out temperature, as well as tF, where F is the res- 
onance width and r is the reinteraction time. Observing 
two ratios where the two particles have the same chemical 
composition, but different m and F, such as yl(1520)/yl vs 
K*/K, or Z'(1385)/yl vs K*/K, could therefore be used 
to extract the magnitude of the freeze-out temperature 
and the re-interaction time. Studies of this type are still 
in progress; As we will show, A{1520) / A and K* /K seem 
to be compatible with sudden freeze-out, provided freeze- 
out happens in a super-cooled over-saturated state [SUHl 
[To] . Other preliminary results, such as p/tt, A/p [8], and 
now S* / A [31j seem to be produced in excess of the statis- 
tical model [51 l30l[5T| . both equilibrium and not. It is diffi- 
cult to see how a long re-interacting phase would produce 
such a result: Resonances whose interaction cross-section 
is small w.r.t. the timescale of collective expansion would 
generally be depleted by the dominance of rescattering 
over regeneration processes at the detectable (on-shell) 
mass range. More strongly interacting resonances would 
be re-thermalized at a smaller, close to thermal freeze- 
out temperature. Both of these scenarios would generally 
result in a suppression, rather than an enhancement, of 
directly detectable resonances. Transport model studies 
done on resonances generally confirm this [5^(55] . Yet the 
only resonance, so far, found to be strongly suppressed 
w.r.t. expectations is the SPS ^1(1520) [34], and even that 
measurement is to date preliminary. 

The observation of an enhanced fJ.'^ fJ-~ continuum around 
the p peak |35] has been pointed to as evidence of p broad- 
ening, which in turn would signify a long hadronic re- 
interaction phase [35] . The absence of a broadening in 
the nominal peak itself, prevents us from considering this 
as the unique interpretation of experimental data. More- 
over, even a conclusive link of broadening with hadronic 
re-interactions would still give no indication to the length 
of the hadronic rescattering period. Nor it would resolve 
the discrepancies pertaining hadronic resonances encoun- 
tered in the previous paragraph; The lack of modification 
in either mass or width, between p-p and Au-Au seen so 
far [30[[3T] is only compatible with the NA60 result pro- 
vided the p is very quick to thermalize, so [30] sees only 
the p s formed close to thermal freeze-out. But, as argued 
in the previous paragraph, that under-estimates the abun- 
dance of the p and other resonances, since the p/w, K* /K, 
and even yl(1520)/yl ratios point to a freeze-out tempera- 
ture significantly above the 100 MeV, commonly assumed 
to be the "thermal freeze-out" temperature in a staged 
freeze-out scenario. 

To help resolve this ambiguity, we aim to directly in- 
fer the magnitude of hadronic reinteraction by combining, 
within the same analysis, an observable sensitive to the 
chemical freeze-out resonance abundance with the direct 



observation of resonances. As shown in [36] . the measure- 
ment of fluctuations of a ratio is such an observable. 

The fluctuation of a ratio N1/N2 can be computed 
from the fluctuation of the denominator and the numera- 
tor[M]K = ((Z\X)2)/(X)): 



'N1/N2 



((Z\A^i)2) {{AN2)') IAN^AN2 



{Ni 



{N2 



- 2- 



(TVi) {N2 



(7) 



Note the appearance of a negative correlation term be- 
tween A'^i and A'2 stemming from a common resonance 
feed-down {A —>■ pi: will be a source of correlation be- 
tween Np and Nt^). 



{AN^AN2) ^Y.^j^i2 {N, 



(8) 



This correlation term is an invaluable phenomenological 
resource, since it is sensitive to resonance abundance at 
chemical freeze-out [3S]; Further reinteraction, provided 
resonances are not kicked out of the detector's acceptance 
region [23], preserve the correlation even when the orig- 
inal resonance stops being reconstructible. Thus, a mea- 
surement of the ratio fluctuation together with the di- 
rectly detectable resonance yield gauges the amount of 
re-interaction between the thermal and chemical freeze- 
out. If the reinteraction period is short, both observables 
should be described by the same statistical model param- 
eters 

The main experimental problem with fluctuation mea- 
surements is the vulnerability to effects resulting from lim- 
ited detector acceptance. This difficulty can be lessened, to 
some extent, by considering "dynamical" fluctuations, ob- 
tained by subtracting a "static" contribution which should 
be purely Poisson in an ideal detector. 



dyn 



(9) 



o'stat, usually obtained through a Mixed event approach 
[37] . includes a baseline Poisson component, which for a 
ratio N1/N2 can be modeled as 



1 



1 



(10) 



as well as a contribution from detector efficiency and kine- 
matic cuts. Provided certain assumptions for the detec- 
tor response function hold (see appendix A of [37]), sub- 
tracting astat from a should yield a "robust" detector- 
independent observable. 

More complicated to deal with are detector acceptance 
effects affecting particle correlations (the probability for 
both resonance decay products to be within the detector 
acceptance region) [23]. These were not corrected for in 
the present study; RHIC data-points do not include this 
term, while at SPS, due to it's large acceptance, this term 
is likely to be less important. 

Related to the time-scale of the thermal freeze-out is 
the question of the total normalization of the system. This 
quantity can be obtain in fits by fitting yields of particles 
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rather than ratios. It is physicaUy important because it 
is connected to the system's volume at chemical freeze- 
out, which is in turn important to gauge the relative im- 
portance of the hadronic phase in the system's dynamics 
[21] . It is also necessary to obtain the thermal energy and 
entropy content of the system [9], and it's scaling with 
centrality. 

On the other hand, normalization introduces an addi- 
tional source of inter-parameter correlation, since it scales 
in the same way as T and 7^. Measuring fluctuations of 
ratios in addition to yields and rations, once again, can be 
used to disentangle these three quantities. Note, from Eqs. 
[7l and [TOl that fluctuations of ratios, while independent of 
the volume fluctuation scale as the inverse of the absolute 
normalization, aj^-^/j^^ ~ {{V) T^)~^. 

Hence an increase in volume makes the yields go up, 
does not affect the average ratios, but decreases the event- 
by-event fluctuation of ratios. An increase in temperature 
generally affects both yields and ratios, increases the yields 
and decreases the fluctuations. An increase in 7^ increases 
yields, baryon to meson ratios and fluctuations. It is there- 
fore not difficult to fix all three parameters with a data- 
sample of relatively few data-points, provided both yields, 
ratios and fluctuations are present. 



2 Fit to RHIC and SPS data 

We have performed fits for both SPS and RHIC energies, 
using publicly available statistical hadronization software 
[381139]. For RHIC, we have used the same data sample as 
in [10], including cr^T/^+ and cr^^"/^_ [H]. For SPS yields, 
we used the same data sample as in [5] , augmented by pre- 
liminary fluctuation measurements of 

and af^^-. + , _^ |42] . It should be noted that the model 
under consideration is far more plausible when applied to 
RHIC than to SPS data, since the fluctuation measure- 
ment in [42] encompasses a large acceptance region with 
non-trivial momentum cuts, that only with the greatest 
hesitation can be imagined as a sub-system in equilibrium 
with a bath. 

Our fit parameters include the normalization (hope- 
fully related to the system "volume" at chemical freeze- 
out), temperature, Xq.sj^ and 7^ ,5. We also require, by 
implementing them as data-points, strangeness ((s ~s) — 
OiO.Ol), charge and baryon number ((Q) / (B) = Z/A± 
0.01) conservation. 

The results of the analysis are shown in table[l]and Fig. 
[TJ An equilibrium analysis (7g = 1, empty blue squares in 
Fig-P) fits the particle yields and ratios but under-predicts 
'^'k/tt many standard deviations. If u'^J^^ is fitted to- 
gether with ratios (red triangles down in Fig. [T]) but no 
yields, it forces the system volume (all other parameters 
being the same within error) to be unrealistically small 
(- 500/m3 at RHIC, - lOOO/m^ at SPS), thereby under- 
predicting particle yields by several standard deviations. 
It is only the addition of 7^ (circles in Fig[T]) that allows 
fluctuations to be driven to a high enough value while 



maintaining sufficiently high volume to describe the par- 
ticle multiplicities, and sufficiently high temperature to 
describe ratios. 

It is important to underline that both yields and fluc- 
tuations contribute to such a precise determination: Equi- 
librium statistical models can describe most yields and ra- 
tios acceptably with 7^ = 7^ = 1, but fail to describe the 
event-by-event fluctuation. Conversely, transport models 
provide an acceptable description of event-by-event fluc- 
tuations [43] , but fail to describe the yield of multi-strange 
particles [44]. 

Unlike what is sometimes asserted, the yl(1520) and 
K* are acceptably described by the statistical model at 
RHIC. The strongest disagreement arises from U* under- 
prcdiction, at the level of 1.5 standard deviations. We 
await for more measurements of resonances such as A and 
p in central collisions before trying to interpret this under- 
prediction. 

The acceptable description of the /1(1520) and K* yield 
using the same freeze-out temperature as the stable par- 
ticles, and the under-prediction of the S* , makes a case 
for the proposition that the re-interaction period between 
hadronization and freeze-out might be not as significant 
as generally thought. However, the current data is not ca- 
pable to rule out such a long-reinteraction period, since 
the crucial fluctuations, of ratios correlated by resonances 
(e.g. vs c^^^'^ ) are still not available at RHIC 

and are not precise enough at for such a falsification at 
SPS. 

It will also be interesting to see if the SPS resonance 
results (preliminary since 2001 [34], shown in the plot but 
not included in the fit) are confirmed: Here, while the K* 
abundance is acceptably described by both the equilib- 
rium and the non-equilibrium model, /1(1520) is consid- 
erably over-predicted by both through the disagreement 
with the equilibrium model is larger, due to the higher 
freeze-out temperature. 

The discrepancies between this fit, and earlier fits with 
Iq IS] are due to data-set choice. The lower temperature, 
and higher 7^ at SPS is due to the E* ,K* and yl(1520) 
resonance yields at RHIC that push for a higher tempera- 
ture, both in equilibrium and non-equilibrium. Until now, 
unfortunately, no published resonance result exists in SPS, 
although /1(1520) and K* results are available in confer- 
ence proceedings These are displayed in Fig. [T] but 
were not used in the fit. 

Qualitatively, the only major difference between SPS 
and RHIC systems is the higher baryo-chemical potential 
(pb), expected due to the higher initial transparency at 
the higher energy RHIC collisions. The introduction of 
7g as a fit parameter makes the freeze-out temperature 
drop to a value compatible with the QGP super-cooling 
hypothesis, and decreases the volume by ~ 30% (a re- 
sult that goes in the right direction to explain the "HBT 
puzzle" [25]). Due to the very different acceptances in the 
experiments at RHIC and SPS, a direct comparison of the 
volumes would not be meaningful. 

The main disagreement between the data and all mod- 
els is the over-prediction of cr"^/" at SPS. Other works [H] 
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Table 1. Best fit parameters at SPS Pb-Pb ^=17.3 GeV and RHIC An- An ^=200 GeV collisions 





RHIC Au-Au Vs=200 CcV 


SPS Pb-Pb Vs=17.3 GeV 




7. /I 


7. = 1 


7. /I 


7, = 1 


T [MeV] 

^lB = 3rin(A,) [MeV] 
Us =T\n{\\-'' [MeV] 
1X13 = Tin A/3 [MeV] 

Is 

Normalization {fm?) 


143.3 ± 2.5 
20.8 ± 3.9 
4.5 ± 1.1 
-0.4 ± 0.1 
1.541 ± 0.001 
1.980 ± 0.158 
1356 ± 395 


159.6 ± 6.9 
23.2 ± 3.4 
5.1 ± 0.8 
-0.7 ± 0.2 
1 

1.083 ± 0.112 
1689 ± 488 


138.2 ± 2.3 
231.2 ± 19.1 
58.6 ± 0.78 
-4.1 ± 1.1 
1.645 ± 0.003 
1.577 ± 0.116 
2540 ± 658 


151.3 ± 3.8 
248.0 ± 25.1 
60.5 ± 9.3 
-8.8 ± 2.1 
1. 

0.815 ± 0.067 
3746 ± 444 
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10 
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10^*5 
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d<N>/dy 



K77t 



p/p 



p/71 — 

^ A/71 



Ratios 




Fig. 1. (color online) Best fit yields, ratios and fiuctuations for RHIC (left panel) and SPS (right panel). Green circles represent 
a fit where 7, was fitted. Red triangles down, 7, = 1 and yields not fitted. Blue squares, 7, = 1 and fluctuation not fitted. 
Red triangles, 7? = 1, normalization determined from the fluctuation and yields not fitted. The magenta diamond refers to a 
calculation with halved yields, as per Eg. 1111 The blue triangle up (right panel only) refers to enhanced A production (iriA = nip) 



have suggested that the low value of cr^J^ is indicative 
of a rich A resonance abundance at chemical freeze-out. 
The preliminary observation of an enhanced A production 
at RHIC [Sj ties in well with this picture. However, this 
study puts this conjecture in doubt: As can be seen from 
the figure, even if one assumes that m/^ = irip, or 80% 
of the photons come from Z\s (blue triangle up in the right 
panel of Fig. [1]) , it is still not enough to account for the 
discrepancy. 

It is intriguing that halving the fluctuation term in Eq. 

m 

1 {{AN,)^) ^ {{AN,)^) JAN^AN,) 
^Ni/N2 = 77 — 7T7T2 1 7T7T2 ^ /t,t \ /t.t \ ■ y^^) 



{Ni 



{N2 



(Ni) {N2 



produces a result compatible with observations (magenta 
diamond in Fig. [T]) for cr^J". This is the scaling expected, 
in the thermodynamic limit, if the ensemble physically rel- 
evant for baryons were canonical rather than GC [T7l[T8] . 
The roughness of this estimation, and the preliminary na- 
ture of the data-points in question prevents us from draw- 
ing conclusions from the result. We also point out that 
if this scaling is applied to Kaons (suggesting C ensemble 
for strangeness) rather than protons (suggesting C ensem- 
ble for Baryon number), the model fails miserably at both 



SPS and RHIC, as Fig. [T] shows. We eagerly await more 
complete and rigorous studies using the canonical ensem- 
ble [19] to further clarify these issues. 

In conclusion, we have used preliminary experimental 
data to show that, at both 200 GeV Au-Au collisions and 
17.6 GeV Pb-Pb collisions, the equilibrium model is un- 
able to describe both yields and fluctuations within the 
same statistical parameters. The non-equilibrium model, 
in contrast, succeeds in describing almost all of the yields 
and fluctuations measured so far at SPS and RHIC, with 
the parameters expected from a scenario where non-equilibrium 
arises through a phase transition from a high entropy 
state, with super-cooling and oversaturation of phase space. 
Some preliminary SPS data-points, such as the yl(152G) 
yield and the c^J", are however significantly over- predicted 
in both equilibrium and non-equilibrium models. We await 
more published data to determine weather the non-equilibrium 
model is really capable of accounting for both yields and 
fluctuations in all light and strange hadrons produced in 
heavy ion collisions. 

Work supported in part by grants from the Natural 
Sciences and Engineering research council of Canada, the 
Fonds Nature et Technologies of Quebec, and the Tom- 
linson foundation. We would like to thank J. Rafelski, C. 
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